Astronomy & Astrophysics manuscript no. terzan6 


February 2, 2008 


(DOI: will be inserted by hand later) 





Bursts, eclipses, dips and a refined position for the luminous 
low-mass X-ray binary in the globular cluster Terzan 6 

J.J.M. in 't F. Hulleman^, C.B. Markwardt^'^, M. Mendez^ E. Kuulkers^, R. Cornelisse^'^, 

J. Heise^'^, T.E. Strohmayer^ , and F. Verbunt'^ 



SRON National Institute for Space Research, Sorbonnelaan 2, NL - 3584 CA Utrecht, the Netherlands 
Astronomical Institute, Utrecht University, P.O. Box 80000, NL - 3508 TA Utrecht, the Netherlands 



m 

^ ; NASA Goddard Space Flight Center, Code 662, Greenbeh, MD 20771, U.S.A. 

^ ! Dept. of Astronomy, University of Maryland, College Park, MD 20742, U.S.A. 

ESA-ESTEC, Science Operations & Data Systems Division, SCI-SDG, Keplerlaan 1, 2201 AZ Noordwijk, the 
Netherlands 



Received, accepted 

(N ■ Abstract. GRS 1747-312 is a bright transient X-ray source in the globular cluster Terzan 6 with quasi-periodic 

^ , outbursts approximately every 4.5 months. We carried out 2-60 keV target-of-opportunity observations during 

eight outbursts with the Proportional Counter Array on the RXTE satellite, for a total exposure time of 301 ks, and 
detect the first unambiguous thermonuclear X-ray bursts from this source. This identifies the compact accretor in 
this binary as a neutron star. The neutron star identification implies that twelve out of thirteen luminous (above 
^— J ' 10^^ erg s~^) X-ray sources in Galactic globular clusters harbor neutron stars, with AC211's nature (in M15) 

I remaining elusive. We observed 24 transitions of eclipses of the X-ray emitting region by the companion star and 

. are able to improve the accuracy of the orbital period by a factor of 10''. The period is P = 0.514980303(7) d. 

We do not detect a period derivative with an upper limit of \P/P\ = 3 x 10~* yr~^. Archival Chandra data were 
analyzed to further refine the X-ray position, and the cluster's center of gravity was re-determined from optical 
data resulting in a correction amounting to 2 core radii. We find that GRS 1747-312 is 0.2 ± 0.2 core radii from 
O ' the cluster center. 
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1. Introduction ray bursts. Such bursts are characterized by fast rises (less 

than a few seconds), slower exponential-like decays with 

GRS 1747-312 is a transient X-ray source that was first e-folding decay times between a few seconds and a few 

detected in 1990 with ROSAT (Predehl et al. 1991) and tens of minutes, and blackbody spectra that exhibit cool- 

Granat (Pavlinsky et al. 1994). Within the uncertainty ing temperatures during the decay (e.g., Lewin et al. 1993; 

of 1' it was found to coincide with the globular cluster Strohmayer & Bildsten 2003). Type-I X-ray bursts are ex- 

(GC) Terzan 6. Later this was confirmed by a more accu- plained as thermonuclear flashes on the surfaces of neutron 

rate ROSAT localization (Verbunt et al. 1995). Terzan 6 stars (Maraschi & Cavahere 1977; Woosley & Taam 1976). 

is a core-collapsed metabrich GC which is offset from the GRS 1747-312 and AC211 in M15 (White & Angelini 

Galactic center by 2?6. The core and half-light radii are 2001) are the remaining cases for which the nature of the 

3'.'3 and 26" (Trager et al. 1995), the distance 9.5^2:5 kpc compact objects needs to be resolved. Knowledge about 

(Barbuy et al. 1997; Kuulkers et al. 2003). It is heavily the nature is of interest for studies about the evolution 

reddened with E{B ~ V) ^ 2.24 (Barbuy et al. 1997) of GCs (i.e., the probabilities of tidal capture and escape 

which is among the highest values of GCs harboring lu- probabilities of the black holes versus the lighter neutron 

minous (> 10^'' erg s^^) X-ray sources. There are thir- stars) , particularly now that there is evidence that GCs in 

teen luminous X-ray sources in twelve globular clusters other galaxies do contain black hole X-ray binaries (e.g., 

in the Galaxy. Thus far, eleven of the thirteen bright GC Angehni et al. 2001 and Sarazin et al. 2001). 

X-ray sources were proven to harbor neutron stars accret- A previous study of mostly moderately sensitive X- 

ing matter from a Roche-lobe filling low-mass companion ray data showed that the outbursts of GRS 1747-312 fol- 

star. The proof was delivered by the detection of type-I X- low each other relatively quickly: In 't Zand et al. (2000) 
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Fig. 1. Time history of GRS 1747-312 photon rate in the 
ASM, between 6 January 1996 and 26 September 2002. A 
bias level of 0.5 c s~^ was subtracted. Sixteen outbursts 
can be recognized. Gaps correspond to times when the sun 
crosses the field of view of the source and observations are 
impossible. The time resolution is 4 d. 



Fig. 2. Time history of GRS 1747-312 photon rate 
during PCA bulge scans (Swank & Markwardt 2001; 
Markwardt & Swank 2002), between 5 February 1999 and 
29 September 2002. Gaps in this plot are also caused by 
sun passages. 



detected five outbursts between 1996 and 1999 with re- 
currence times of about half a year. It was also estab- 
lished that GRS 1747-312 is one of now six known low- 
mass X-ray binaries that exhibit complete X-ray eclipses; 
the others are EXO 0748-676 (Parmar et al. 1986, Wolff 
et al. 2002), MXB 1658-298 (Cominsky & Wood 1984, 
Oosterbroek et al. 2001), Her X-1 (Tananbaum et al. 
1972), XTE J1710-281 (Markwardt et al. 1999) and AX 
J1745.6-2901 (Maeda et al. 1996; Kennea & Skinner 1996). 
The measured orbital period for GRS 1747-312 is 12.360± 
0.009 hr, the eclipse duration 43 min. 

We initiated a program of dedicated and sensi- 
tive target-of-opportunity (ToO) X-ray observations with 
RXTE to search for signatures regarding the nature of 
the accretor in GRS 1747-312, in particular type-I X-ray 
bursts which can unambiguously be identified as ther- 
monuclear flashes in a neutron star atmosphere (like for 
the other eleven previously mentioned GC cases), and to 
measure more eclipses to improve the knowledge about the 
orbital period. The present paper concerns the full analysis 
of all X-ray bursts, the timing study of the eclipses, and a 
discussion of the recurrence behavior of the transient out- 
bursts. Furthermore, we analyzed archival Chandra data 
to determine the most accurate X-ray position thus far. In 
connection to this we re-determined the cluster's center of 
gravity from optical data and find that the X-ray source 
is within one core radius from the center. 

2. ASM and PCA long-term light curves 

The All-Sky Monitor (ASM) on RXTE monitors GRS 
1747-312 for about ten months every year, for up to a few 
tens of times per day during 90 sec dwells. The sensitivity 
is typically a few tens of mCrabs"'^ per day of observations 
(Levine et al. 1996). Although limited, the sensitivity is 
sufficient to pick up outbursts by GRS 1747-312 which 
typically peak at 30 mCrab. The resulting light curve is 
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Fig. 3. Same as Fig.|21 except that the times of the 8 max- 
ima have been subtracted (different symbols for different 
outbursts). The vertical lines indicate the times of detec- 
tions of type-I X-ray bursts with the PCA (solid) and of 
burst- like features with the WFCs (dashed) . The three low 
points during the outburst phase have been verified to be 
caused by eclipses. 



presented in Fig. ^ Sixteen outbursts are readily recog- 
nized by eye, a few (maybe two) were probably missed at 
times when the sun passed by GRS 1747-312 and obser- 
vations were impossible. 

Since February 1999 the Proportional Counter Array 
(PCA, see next section) is employed in slow semi-weekly 
scans of a 16° x 16° field centered on the Galactic center 
(Swank & Markwardt 2001). Fluxes of all point sources 
in this field are monitored, including that of GRS 1747- 
312. The observation frequency is one to two orders of 
magnitude less than that of the ASM observations, but 
this is compensated by a considerable gain in sensitivity. 
Fig. 121 shows the resulting light curve. 

If we assume that two outbursts were missed in the 
ASM data (between outbursts VIII and IX and between 
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XII and XIII), the average recurrence time is 136 days. 
However, the recurrence time varies between 76 days (be- 
tween outbursts VII and VIII, see Fig. ^) and 171 days 
(between outbursts XI and XII). There is a clustering of 
recurrence times in the 130-142 day range: 7 out of 13 
securely determined cases are in this range. 

Nearly all outbursts have similar shapes and peak in- 
tensities of about 380 c s'^ per 5 PCUs or 2 ASM c s'^, 
corresponding to 0.03 Crab units; only two outburst are 
significantly different: the peak intensity of outburst VII 
is half the usual value and outburst XII is 25% larger. 
These deviations are accompanied with extremes in the 
recurrence time. Figure |3| shows the PC A scan light curve 
folded with respect to the times of maximum of the seven 
outbursts. All PCA-detected outbursts last 4 weeks except 
for the last outburst which had low-level activity preceed- 
ing the maximum for one week. 

The average intensity from outburst to outburst may 
accurately be determined from those outbursts that were 
well sampled by the PCA bulge scans and were accompa- 
nied by unambiguous measurements of the recurrence time 
since the previous outburst. Five outbursts meet these cri- 
teria: X, XI, XII, XV and XVI (including the deviating 
outburst XII) . We determined the photon fluence per out- 
burst in a model-independent manner by 1) subtracting a 
background level as determined from 3 points preceed- 
ing the peak by at least 7 days; 2) multiplying each data 
point with the interval time to the next data point; and 
3) summing all values. The resulting fluences divided by 
the recurrence time provides average intensities between 
outbursts ranging from 50± 2 to 58± 1 c s~^ per 5 PCUs. 
The fractional rms value over the 5 measurements is 5%. 



3. RXTE-PCA pointed observations 

The Proportional Counter Array (PCA; for a detailed de- 
scription, see Jahoda et al. 1996) on RXTE consists of an 
array of 5 co-aligned Proportional Counter Units (PCUs) 
that are sensitive to photons of energy 2 to 60 keV with a 
total collecting area of 6500 cm^. The spectral resolution 
is 18% full- width at half maximum (FWHM) at 6 keV and 
the field of view is 1° FWHM. During any observation the 
number of active PCUs may vary between 1 and 5. 

The data on GRS 1747-312 were obtained in 80 ToO 
observations that cover eight outbursts in 1998-2002, see 
Tabled The total exposure time is 301 ks. Thirteen obser- 
vations (OBSIds 50034-01-05 through 50034-01-06) were 
scheduled to search for an anticipated outburst, because 
the usual trigger observations by the ASM and PCA bulge 
scans were absent due to small angular distances to the 
Sun. Unfortunately, that outburst did not materialize. On 
four outbursts wide coverage was obtained (Feb 2000, Jul 
2000, June 2001 and October 2001). Coverage of the other 
four outbursts was limited to one or two observations of 
eclipses. Many of the individual exposures were scheduled 
to cover either an eclipse ingress or egress and not the full 
eclipse. An eclipse is rather long and spending much of 



Table 1. Summary of pointed RXTE observations. The 
dividing lines indicate different outbursts. 



Out- 


ObsID 


Dates 


No. 


Exp. 


burst 






Obs. 


(s) 


VIII 


30427-01-01 


1998-09-13 


1 


3283 


IX 


40419-01-01 


1999-05-12 


1 


3356 


X 


40419-01-02 


1999-09-27 


2 


6592 




40419-01-03 


1999-10-13 


1 


3584 




40419-01-04 


1999-10-20 


1 


2432 


XI 


50034-01-01 


2000-02-18,20,21,24 


7 


26546 




50034-01-02 


2000-02-25,26,27,28 


9 


34272 




50034-01-03 


2000-03-11,12 


2 


2908 




50034-01-04 


2000-03-15 


1 


3518 




50034-01-05 


2000-06-30 


7 


11886 






2000-07-1,2,3,4,5,6 








50034-01-06 


2000-07-7,8,9,10,11 


6 


7791 






2000-07-12 






XII 


50034-01-07 


2000-08-10 


2 


13521 




50034-01-08 


2000-08-11,12,15 


3 


5717 




50034-01-09 


2000-09-09,11 


3 


6000''^ 


XIII 


60038-01-01 


2001-05-31 


1 


13860 




60038-01-03 


2001-06-03,06 


2 


5609 




60038-01-04 


2001-06-10,12,13,16 


4 


20942 




60038-01-05 


2001-06-17 


2 


9300 




60038-01-06 


2001-07-1,2,5 


3 


6556 


XIV 


60038-01-02 


2001-10-21,23,24 


8 


43983 




60038-01-07 


2001-10-26,27,29,30 


8 


48682 






2001-10-31 








60038-01-08 


2001-11-02,6,8 


3 


14307 




60038-01-09 


2001-11-09 


1 


6182 


XV 


60038-01-10 


2002-03-14 


1 


1377 


XVI 


60038-01-11 


2002-07-16 


1 


1101 


t These 


observations 


had bad attitude 


data 


and were 



disregarded. 



the exposure budget on that would have diminished the 
chances to detect an X-ray burst. 

4. X-ray bursts 

4.1. Detections 

In outbursts XIII and XIV we were succesful in catch- 
ing bursts. They were detected on June 17 (burst no. 1), 
October 29 (no. 2), November 2 (no. 3) and 6, 2001 (no. 
4) ; all occurred during the decay phase of an outburst (see 
Fig.|2Il. While the bursts all have durations of 10 to 15 sec, 
they rather vary in profile and peak flux (see top panels 
of Fig. ^ and there is a slight concern that other sources 
in the field of view are responsible. Therefore, we applied 
a localization procedure to those bursts that employs the 
slight offsets that exist between PCUs (Strohmayer et al. 
1997). The result is consistent, at the 90% confidence level, 
with all bursts coming from within 0?5 of the position of 
GRS 1747-312. No other bursters are known within this re- 
gion, although SAX J1752. 3-3138 (Cocchi et al. 2001) is at 
the edge. Furthermore, no otherwise (transiently) bright 
X-ray sources are known in this region. Given that the 
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Fig. 4. Time-resolved spectroscopy of the bursts. The upper panel of each burst shows the raw full-bandpass 
background-subtracted (including the persistent emission from GRS 1747-312) photon count rate, the 2nd panel from 
above the bolometric absorption-corrected flux in units of 10~^ erg s~^cm~^, the 3rd panel the blackbody temperature 
in units of keV, the 4th panel the radius of the spherical blackbody emission area in km (for a distance of 9.5 kpc), and 
the bottom panel the reduced values for the spectral fit to each time-resolved spectrum. We note that the orbital 
phases with respect to mid-eclipse for these bursts are 0.86, 0.28, 0.21 and 0.49 (see Sect. IS)). The abscissa reference 
times are in UTC. 



bursts all occurred when GRS 1747-312 was in an active 
state, we conclude that they are from the same source. 



Between mid 1996 and mid 2002 the Wide Field 
Cameras (WFCs; Jager et al. 1997) on BeppoSAX (Boella 
et al. 1997) were used to monitor a 40° x 40° field at 5' 
resolution around the galactic center during Feb- Apr and 
Aug-Oct of each year, with exposure times of 100 ksec to 
1 Msec and observation frequencies of 4 to 67 per visibil- 
ity window (e.g.. In 't Zand 2001). The total net exposure 
time is 7 Msec. The bandpass of the WFC is 2 to 28 keV, 
and the typical sensitivity is 10 mCrab in one observation. 
Part of these data were the basis of the detection of the 
12.4 hr orbital period of GRS 1747-312 (In 't Zand et al. 
2000). The WFC observations revealed the detection of 
three burst-like features, with durations of less than 10 s 
and peak fluxes comparable to those of the brightest PCA- 
detected bursts. The bursts occurred during the decays of 
outbursts XI and XII (see Figs. |21 and |3J); two were within 
half a day from each other. Although the data are of in- 
sufficient statistical quality to allow a detailed spectral 
analysis, they do confirm with high localization precision 
the bursting nature of GRS 1747-312. 



4.2. Time-resolved spectroscopy 

We carried out a time-resolved spectral analysis of the 
four PCA-detected bursts. Preliminary results were pre- 
sented in Kuulkers et al. (2003). Event mode data with 
a resolution of 2^^^ sec and 64 energy channels were 
binned to 0.25 through 4 sec, the background was sub- 
tracted as well as the pre-burst persistent emission from 
GRS 1747-312, a correction for deadtime was applied, and 
the spectral data were, between 2 and 20 keV, fitted with 
a single-temperature blackbody model and constant ab- 
sorption {Nn = 1.4 X 10^^ cm~^, see In 't Zand et al. 
2000). Furthermore, spectral channels were combined so 
that each spectral bin contained at least 15 photons, to 
ensure Gaussian statistics and the applicability of the 
statistic. The smallest number of spectral bins was 3. The 
results are summarized in Fig. ^ together with the raw 
photon count rate profiles. The spectral fits are generally 
acceptable. All bursts show a temperature profile that de- 
creases during the decaying phase. The blackbody spec- 
trum and the cooling during decay uniquely identify the 
bursts as type-I X-ray bursts which are explained as due 
to thermonuclear flashes in the upper layers of neutron 
stars (see §1). 

Two of the bursts show a time profile for the emis- 
sion area which is not constant and anti-correlated with 
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temperature variations. This is indicative of photospheric 
radius expansion (by a factor of 3 for burst 2 and 6 for 
burst 3) due to near-Eddington fluxes during these events. 
Remarkably, the unabsorbed bolometric peak flux during 
these two bursts differs by a factor of 1.8. 

4.3. Search for burst oscillations 

We searched the four X-ray bursts for osciflations. Light 
curves were produced at the full time resolution of the 
data, 2~^^ s, starting ~ 5 s before the onset of the bursts 
and extending up to times when the intensity had decayed 
back to within ~ 10% of the preburst value. We used data 
for the whole energy band, nominally 2 — 60 keV, and from 
all the detectors that were on during each burst (three de- 
tectors during burst 1, 3, and 4; and four detectors during 
burst 2). We calculated one Fourier power spectrum of the 
whole burst for each burst, and power spectra of 1, 2, and 
4-s intervals covering the whole duration of the bursts. In 
the case of the 1,2, and 4-s power spectra, to increase the 
sensitivity to transient oscillations lasting less than the 
time length of the power spectra, we set the beginning of 
the interval for which we calculated the power spectrum 
0.125, 0.250, and 0.5 s, respectively, after the beginning 
of the previous interval (the power spectra were therefore 
not independent). For the brightest burst (no. 3), we also 
calculated power spectra of 1) 0.25 and 0.5 s intervals, 
shifting by 1/32 and 1/16 s at a time; and 2) for > 6 keV 
photons only. 

We did not flnd any significant oscillation in the 
0.025 — 1024 Hz frequency range in any of the bursts. 
Because of the rapid change in intensity during the bursts, 
we calculated upper limits separately for intervals during 
which the source intensity was larger than 90% of the 
peak intensity, and the rest of the data. For the faintest 
burst (no. 1) the observed count rate was a factor of ~ 7 
lower than for the brightest burst (no. 3), and therefore 
the upper limits vary considerably from burst to burst. 
The 95 % confidence upper limits for oscillations (for the 
method, see Vaughan et al. 1994) in bursts 1 and 4 dur- 
ing the bright and weak parts are 29.8% rms and 63.8% 
rms, respectively. For bursts 2 and 3, the 95% confidence 
upper limits in the bright and weak parts are 11.5% rms 
and 47.3% rms, respectively. 

In the 11 sources that currently are known to exhibit 
burst oscillations, the oscillations may reach high ampli- 
tudes (>50% rms) during short intervals (typically 0.1 s) 
if they occur during the rising phase and low amplitudes 
(<15%) during the decaying phase (for a recent review, 
see Strohmayer & Bildsten 2003). Therefore, our upper 
limits are not very constraining. 

5. Eclipses and dips 

Twelve eclipse ingresses and twelve egresses were ob- 
served, in twenty-one eclipses distributed over seven out- 
bursts. Two eclipses were observed completely; for another 
one the ingress and egress was covered with a data gap 
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Fig. 6. Observed minus calculated ('0-C') diagram for the 
eclipse ingresses (circular symbols connected by a dotted 
line) and egresses (rectangular symbols connected by a 
dashed line) as observed with RXTE, with respect to the 
respective ephemerises (see text). 

during the eclipse; for the rest only the ingress or egress 
was covered. Timing of the transitions provides the means 
to improve the accuracy of the orbital ephemeris signifi- 
cantly. The orbital period thus far had only a 31 sec ac- 
curacy (In 't Zand et al. 2000). More accuracy enables 
studies of changes in the orbital period. To increase the 
time baseline, we include the timing of the egress observed 
with the BeppoSAX Narrow Field Instruments on 1998 
September 6 (In 't Zand et al. 2000). Thus, the time base- 
line is 3.9 yr. 

Ingress and egress last between 30 and 40 sec. We took 
the mid points between minimum and maximum flux as 
a reference for the timing of the eclipse transitions. The 
transition profiles appear too variable to trust methods 
that involve profile fitting. Therefore, we resorted to the 
following procedure. As a baseline, we used standard 1 
data which provides 0.125 s resolution. The times were 
corrected to the solar system barycenter^. We identified 
the eclipses in a light curve of the whole data set at a res- 
olution of 1 s and selected data stretches of 200 s centered 
on the initial midpoint determinations (by eye) of ingress 
or egress. The light curves of the center 120 s parts of these 
data stretches are displayed in Fig. |S1 One ingress had to 
be excluded from the analysis (at MJD 52060) because the 
data stream stopped before total eclipse was reached; one 
egress was excluded because the source was too faint (at 
MJD 51601). We determined the weighted average of the 
initial and final 50 s of each data stretch, except for the 

using the f axbary routine in ftools version 5.2 
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Fig. 5. Light curves of 120 s data stretches around 24 echpse transitions in GRS 1747-312. The abscissa refers to time 
in second, with a reference time that is indicated above each plot in MJD (using Terrestrial Time and corrected to 
solar system barycenter; TDB), and the ordinate shows the photon rate normalized for 5 PCUs (note that the actual 
number of active PCUs is not constant over all 24 data stretches and is indicated at the top of each plot). 



very last egress in which a small flare-like feature needed 
to be avoided (see Fig. EJ- This excludes the eclipse tran- 
sitions. The time when the flux reaches the halfway point 
between these two averages is determined by identifying 
the data point closest to this halfway flux, selecting the 17 
data points nearest to this point, making a least-squares 
linear fit to these 17 data points, and finally calculating 
the interpolated time of the halfway point. This procedure 
was developed more or less empirically. It avoids compli- 
cations due to changing transition profiles, particularly 



near the start and end of each transition. There may be a 
bias introduced by the variability of the source outside the 
eclipse but we estimate this to be negligible compared to 
the statistical accuracy. The statistical errors are of order 
1 s (la) in this analysis. The results of the timing are given 
in Table 121 In addition to the timing of the mid-transition 
points, we determined a measure for the transition dura- 
tion by calculating the time it takes the linear fit to bridge 
the weight averages of the initial and final 50 s stretches. 
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Table 2. Times and durations of all eclipse transitions 
observed thus far with sufficient sensitivity 



-30 -20 -10 10 20 
Time since mid transition [s] 

Fig. 7. Average ingress (thick) and egress profile (thin) 
after re-normalization of both and mirroring the time axis 
of the egress profile. The average profiles have been con- 
structed by aligning individual profiles at the time of mid 
transition (see Table 0). 

From these data we determine the ephemeris for the 
ingresses and egresses. The deviations from these solu- 
tions (the "O — C" values) are presented in Table |21 
and Fig. El Ten out of the twenty values are larger 
(in a positive or negative sense) than 3ct, but there 
are no clear trends in O — C and we conclude that 
there is no evidence for systematic changes in the or- 
bital period. Rather we suspect that the larger-than- 
statistical errors are due to imperfections of the method. 
Assuming this is the case, we multiply the errors of 
the linear solutions of the ingresses and egresses by the 
square-root of the Xicd values {\/K6 for the ingresses 
and Vis. 7 for the egresses), thereby effectively forcing 
Xrcd to 1. The ephemeris for the ingresses is Ti(n) = 
MJD (TDB) 52066.259473(5) + 0.514980311(8)n and for 
the egresses T^{n) = MJD (TDB) 52066.289497(10) + 
0.514980275(15)n. The difference in the two derived or- 
bital periods is 3.1 ms with an error of 1.7 ms. We take for 
the orbital period the weighted mean of 0.514980303(7) d. 
This is only 0.15 s shorter than the previous determi- 
nation (In 't Zand et al. 2000). We determine an upper 
limit to the period derivative of \P/P\ = 1 x 10^^^ s^^ or 
3 X 10^* yr~^. The average eclipse duration is 2596(2) s 
(or 43.27 min, or 0.05834 orbital periods). 

It may be that the larger-than-statistical errors in the 
mid-transition times are due to intrinsic properties rather 
than imperfections of the method. A similar phenomenon 
has been observed in EXO 0748-676. In addition to a long 
term period evolution trend, a significant amount of jitter 
is detected in that source (Wolff et al. 2002). Wolff et al. 
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-1.2 (-2.2) 


33.6 + 3.4 


52066.259477 


0.47 





+0.3 (+0.7) 


29.7 + 2.6 


52077.074080 


0.45 


21 


+1.8 (+3.9) 


25.5 + 2.3 


52077.589014 


1.42 


22 


-2.2 (-1.6) 


31.3 + 5.2 


52203.244242 


0.36 


266 


+0.5 (+1.5) 


25.9 ± 1.8 


52347.438718 


0.47 


546 


-0.4 (-0.9) 


21.6 + 2.0 


egresses 


51062.592982'' 


3 


-1949 


+3.6 (+1.2cr) 




51594.567535 


0.52 


-916 


-2.6 (-5.0) 


31.6 + 3.2 


51600.747403 


1.14 


-904 


+6.4 (+5.7) 


39.3 + 8.3 


51602.807303 


4.08 


-900 


+4.6 (+1.1) 


68.1 + 22.1 


51615.166779 


4.68 


-876 


+0.2 (+0.1) 


50.5 + 10.0 


51768.115903 


0.50 


-579 


-1.3 (-2.6) 


36.1 + 3.6 


52060.624742 


0.60 


-11 


+2.4 (+3.9) 


36.2 + 4.1 


52066.289531 


0.95 





+2.9 (+3.1) 


46.9 + 6.5 


52073.499250 


0.97 


14 
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t Number of orbits: 
text for definition; 



* Deviation with respect to linear fit; 
" BeppoS AX-NFI observation 



demonstrate that the eclipse phasing of EXO 0748-676 is 
compatible with a random walk process. They speculate 
that the companion star may become significantly tidally 
distorted, enough so that the times of eclipse ingress and 
egress can be shifted from their nominal value. The long 
term period evolution may be driven by spin-orbit cou- 
pling via the quadrupole moment of the companion star. 
In the present analysis we are only interested in the long- 
term orbital period evolution of GRS 1747-312 and we 
regard the jittering as noise. 

The average transition profiles are given in Fig. |7| 
Surprisingly, there is a difference between the average 
ingress and egress. In particular, the average egress ap- 
pears to progress about 1.3 times slower than the ingress. 
This is confirmed by the averages of the durations as given 
per transition in Tabled they are 27.3 ± 0.8 s for ingress 
and 33.8 ± 1.3 s for egress. 

In Fig.|Slwe present the average ingress and egress pro- 
files for two separate bands: 2-5 keV and 5-20 keV. These 
show that the eclipse lasts approximately 25 s longer in 
the low-energy band, indicating absorption in the atmo- 
sphere of the obscuring companion star. This is spectrally 
confirmed by the BeppoSAX measurement (In 't Zand et 
al. 2000). 

We tested the transition profiles in the two band- 
passes with a simple model involving an isothermal spher- 
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Fig. 8. In thick lines, the average ingress profile (left) and 
egress profile (right) in the 2-5 keV band; in thin lines in 
the 5-20 keV band. The averages have been taken over 
6 good-quality ingresses and 8 egresses. Background sub- 
traction consisted of two steps: first for particle-induced 
and cosmic background (using the pcabackest tool in 
f tools version 5.2) and second for residual emission dur- 
ing total eclipse. The 5-20 keV intensities have been renor- 
malized to match the plateau values of the 2-5 keV curves. 
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ical atmosphere in hydrostatic equilibrium where the den- 
sity is an exponential function of height, assuming a con- 
stant source spectrum that is equal to that measured with 
BeppoSAX (In 't Zand et al. 2000). This model turns out 
to give a good representation of the data. Both the low 
and high-energy profiles are well reproduced as is the lag 
between them. We find neutron star traverse times over 
the scale height of 6.0 ± 0.5 and 8.5 ± 0.5 s for the aver- 
age ingress and egress respectively, with goodness-of-fits 
of x?ed = 0-893 (400 dof) and 0.989 (386 dof) respectively. 
With the eclipse duration this results in upper limits for 
the true scale height in terms of stellar radius of 0.2 to 
0.3%. We also tested hnearly extended sources instead of 
a point source, as would be expected for a nearly edge-on 
accretion disk, and find fits that are only marginally better 
by about 3% in Xredi with source-size equivalent traverse 
times of 14 and 20 s and scale heights only 0.5 to 1.0 s 
smaller than for a point source. We regard the evidence 
for an extended source as weak. 

As in other eclipsing low-mass X-ray binaries (LMXBs; 
EXO 0748-676, MXB 1658-298, XTE J1710-281), GRS 
1747-312 shows dipping activity (for an example, see 
Fig.O. Dips are thought to be caused by coverage of the 
emission region by irregular thickening of the accretion 
disk (White & Swank 1982). Thanks to the large inclina- 
tion angle, inferable from the occurrence of eclipses, the 
presence of dips is not unexpected. Dipping activity was 
observed during outbursts XI (during two different binary 
orbits), XII (one orbit) and XIV (one orbit). The orbital 
phase with respect to mid-eclipse was between 0.58 and 
0.74. The dip depth never exceeds 70%, indicating that 
the emission region is never completely obscured by the 
disk irregularities. The dips last typically a fraction of a 
minute. 



Fig. 9. Upper panel: 2-20 keV light curve during a period 
of dipping activity; lower panel: 5-10 to 2-5 keV hardness 
ratio. The times are solar-system barycentered and the 
photon rates have been background subtracted. The back- 
ground count rate includes particle-induced background, 
diffuse background and the local Galactic background as 
determined from the full eclipse observed on MJD 52206 
at which time the same PCUs were employed. 



6. Accurate positions of X-ray source and cluster 
center 

Chandra observed GRS 1747-312 on March 9, 2000, from 
5:44:44 to 8:41:16 UT (PI J. Grindlay). This observa- 
tion was triggered by our RXTE observations showing the 
source to be active. The High Resolution Camera HRC-I 
was used for a total net exposure time of 9964 s. No grat- 
ing was used. The HRC-I is a microchannel plate detector 
with a field of view of 30'x30' and an angular resolu- 
tion of 0'.'4 (FWHM), a bandpass of 0.1 to 10.0 keV with 
227 cm^ effective area at 1 keV, and minimal spectral res- 
olution (Murray et al. 2000). We extracted the data from 
the public archive for an analysis of the position. Such 
an analysis was performed previously by Revnitsev et al. 
(2002). However, we anticipated to improve on that by 
attempting an astrometric solution of the image through 
X-ray detections of cataloged stars. 

GRS 1747-312 was significantly detected. A total of 
15x10"^ photons were detected and the intensity was 
1.50±0.01 s^^. Slow variability was observed with a « 20% 
modulation around the average on a time scale of 10'^ s. 
Applying the 'celldetect' tool from CIAO (version 2.2) to 
the image, we find four more point sources with a net num- 
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Fig. 10. Gunn z plate of Terzan 6 (Barbuy et al. 1997). 
The weakest stars are approximately z = 23. The position 
of the bright transient is marked by a white circle (with 
error radius of O'.'S), the newly determined position of the 
center of gravity is indicated by two white lines, and the 
old position by a black cross. 



ber of photons in excess of 10, and three sources between 
6 and 10. We determine the pixel positions of all eight 
sources. Using the satellite aspect solution and requiring 
a match within 1" (the absolute astrometric accuracy is 
quoted by the Chandra team as O'.'S to O'.'S at 90% con- 
fidence), we tentatively identify three X-ray detections of 
stars in the USNO-Bl.O catalog (Monet et al. 2003). They 
are USNO0587-0575246 (&usno = 12.46, rusNO = 11.45), 
USNO0587-0575103 (14.79 and 13.32) and USNO0586- 
0575381 (16.10 and 14.19). From the same optical im- 
age we determined that the probability for chance coin- 
cidence of a star brighter than rusNO = 14.19 with one 
l"-radius circle is 7 x 10^* and we conclude that the iden- 
tifications are correct. We determined (X-ray minus op- 
tical) offset corrections of 0'.'68 ± 0'.'12 and 0'.'53 ± 0'.'12 
in a2ooo.o and (52ooo.o respectively. Applying these cor- 
rections, we arrive at a position for GRS 1747-312 of 
^2000.0 = 17'^50'"46'.'862, (^2000.0 = - 31°16'28'.'86 with 
an error radius of 0'.'4 (95% confidence). This position is 
consistent with the position determined by Revnitsev et 
al. (2002), but has an 1.8 times better accuracy (when 
compared at the same confidence level). 

On 17 May 1994, optical images of Terzan 6 were ac- 
quired with the 3.55 m New Technology Telescope (NTT) 
of the European Southern Observatory, using the SUST 
1 camera which employed a 1024x1024 CCD at the 
Nasmyth focus B, giving a plate scale of 0'.'13 per pixel 
(Barbuy et al. 1997). These images were taken under ex- 
cellent seeing conditions (0'.'6) in / (300 s), V (480 s) and 
Gunn z (120 s) filters. We carried out astrometry of the V 



and z plates relative to the USNO-Bl.O catalog. Centroids 
of observed star profiles were determined by Gaussian fit- 
ting. The astrometric solution consisted of a zero point 
position, plate scale, and position angle on the sky. 14 
catalog stars were identified in the V image with devi- 
ations less than 0'.'4, and 19 stars in the z image with 
deviations less than O'.'B. The rms values of the measured 
star positions with respect to the cataloged positions is 
0'.'20 in right ascension and 0'.'14 in declination (0'.'25 and 
0'.'24 respectively for z). Both solutions are less than 0'.'15 
apart. 

The position of the center of Terzan 6 that is used 
throughout the literature (a20oo.o = 17^50™46'.'44, 
<520oo.o = — 31°16'31'.'4) is due to the original discov- 
ery plate (Terzan 1968, see also Terzan 1971). From com- 
parison with Fig. ^] it is clear that this position is sub- 
stantially off the densest concentration. Therefore, we re- 
determined the center through the iterative centroiding 
method introduced by Auriere (1982; see also Picard & 
Johnston 1994). In this method the (unweighted) average 
position is taken over all resolved stars within a certain ra- 
dius from a trial position. This average becomes the trial 
position in the next iteration. The procedure is repeated 
until the average converges. For Terzan 6, a radius of 15" 
was chosen on the basis of the cluster's appearance in our 
image. The resulting center is 02000. = 17'^50™46'.'854, 
(52000.0 = - 31°16'29.38'.'4 with an uncertainty of 0'.'6 as 
evaluated by testing different values for the radius. The ul- 
timate determination of the center is based on the average 
position of 219 stars. This newly determined center is 5'.'7 
away from the old position. We note that this correction 
does not affect the structural parameters of the cluster 
because those were determined completely from relative 
positions (Trager et al. 1995). 

The Chandra position of GRS 1747-312 is 0'.'5 from 
our determination of the center which is equivalent to 0.2 
core radii. This confirms the suggestion from the sub-day 
orbital period and the moderately high collision number 
(Verbunt 2001) that this binary was formed due to tidal 
capture or exchange encounter. 

No sources belonging to Terzan 6 other than GRS 
1747-312 were unambiguously detected in the Chandra 
image. Outside the point spread function of GRS 1747- 
312, the upper limit for the intrinsic luminosity for a 
0.3 keV black body spectrum with A^h = 1-4 x 10^^ cm~^ 
is approximately 1.3 x 10"^'^ erg s~^. The point-spread func- 
tion is roughly given by a Gaussian profile with a standard 
deviation of 0'.'4. It degrades the sensitivity to faint sources 
by a factor of more than two for distances closer to GRS 
1747-312 than 1'.'8. 

7. Discussion 

7.1. Outburst behavior 

The 2-10 keV peak fiux is mostly about 
7xl0~^° erg cm~^s~^. Broad-band measurements 
during outburst VIII allowed an extrapolation to 0.1- 
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200 keV of 1.0x10-9 erg cm.-^s-^ (In 't Zand et al. 
2000). The outburst duration is always 4 weeks. The 
fractional rms in the average intensity from outburst to 
outburst is 5%. If we assume that the spectrum does 
not change drastically from outburst to outburst, the 
luminosity averaged over all times (outburst as well as 
quiescence) is l.SxlO'^^ erg s"^ and the equivalent mass 
accretion rate, assuming 100% efficiency of transform- 
ing the liberated gravitational energy to radiation, is 
about 1.2x10-10 Mgyr"! (i.e., 1% of Eddington). One 
accurately measured outburst (XII) was delayed by 25% 
of the average wait time. The peak of that outburst 
was also 25% larger while the duration was the same, 
so that the average luminosity remained constant. It 
must be a rather constant mass transfer rate from the 
companion star to the accretion disk which is responsible 
for the regular outburst behavior. A similar behavior 
is observed in the Rapid Burster in the GC Liller 1: it 
shows a quasi-periodic outbursts whose period of about 
210 days suddenly switched to 100 days in Nov. 1999 in 
combination with a factor-of-two decrease in the peak 
flux (Masetti 2002). Another well-documented case is 
Aql X-1, with at least 31 outbursts detected thus far (e.g., 
Simon 2002) , but its recurrence time is not so well defined 
and there does not seem to be as strong a correspondence 
between recurrence time and fluence per outburst. 

The transient nature that many LMXBs exhibit is gen- 
erally thought to be related to thermal-viscous instabili- 
ties in the accretion disk where the thermal conditions 
are predominantly determined by heating through X-rays 
from the central object (e.g., Smak 1983; Osaki 1996; van 
Paradijs 1996; Lasota 2001, Simon 2002). These instabil- 
ities become important and the source becomes transient 
if the average X-ray luminosity is below a critical value 
(van Paradijs 1996, see also White et al. 1984). The on- 
set of an outburst occurs if the surface column density 
at any location in the accretion disk exceeds a critical 
value. As a consequence the temperature rises, the vis- 
cosity increases, and heating fronts travel in and outward 
through the disk. Depending on the radius of the ignition, 
differently sized outbursts may result (outward fronts are 
expected to be less efficient, possibly resulting in smaller 
outbursts if the igition occurs in the inner parts of the 
disk). Furthermore, the ignited parts of the disk may ir- 
radiate other cool parts of the disk and ignite those. The 
increased viscosity results in a depletion of the accretion 
disk which may be complete or partial. The latter is more 
likely to happen if the ignition radius is small. The recur- 
rence time, peak flux and outburst duration depend on 
quite a number of parameters and effects, but it seems 
obvious that the quasi-periodicity at least testifies to a 
rather stable mass transfer rate from the secondary. The 
similar recurrence times of GRS 1747-312, Aql X-1 and 
the Rapid Burster probably reflect similar orbital peri- 
ods (19 hrs for Aql X-1, unknown for the Rapid Burster), 
similar mass transfer rates and the fact that they are all 
neutron star systems without considerable truncation of 



the disks that may result in decades-long recurrence times 
such as in black hole systems (e.g., Dubus et al. 2001). 

7.2. Eclipse timing and evolutionary constraints 

As illustrated by Horne et al. (1985), the eclipse duration 
expressed as a fraction of the orbital period constrains 
the mass ratio and inclination. For the low-mass X-ray 
binary in Terzan 6 we find that the inclination must be 
higher than 74? 5, and the mass of the mass donor must 
be higher than 0.07 times the mass of the neutron star. 
Thus, for a canonical 1.4 M0 neutron star, the lower limit 
on the mass of the secondary is 0.1 M0; a likely upper 
limit is 0.8 Mq. Kepler's third law yields a semi-major 
axis of 3.0 to 3.6 Kq for these masses. This implies an 
average orbital velocity of the secondary with respect to 
the primary between 296 and 356 km s^^. 

The eclipse timing resolves an orbital period that has 
remained constant within \P/P\ < 3 x 10~* yc^^ for al- 
most 4 years. This corresponds to a lower limit to the time 
scale of any orbital period change of 3 x lO'' yr which is 
mildly constraining. Only one LMXB has evidence of a 
shorter time scale: P/P = 2 x lO"^ yr in EXO 0748-676 
(Wolff et al. 2002). 

We may compare; these observational constraints with 
theory as follows. In a globular cluster, the most massive 
star still on the main sequence has a mass of about 0.8 Mq. 
A main-sequence star of this mass, or lower, cannot fill its 
Roche lobe in an orbit of 0.515 d. Thus, the donor of GRS 
1747-312 must be larger than a main-sequence star of the 
same mass. If it is larger due to stellar evolution into a 
(sub)giant, then its initial mass was the turnoff mass of 
0.8 Mq. For a neutron-star mass of 1.4 Mq and a compan- 
ion of 0.8 M0, the orbital period gives a distance between 
the two stars of 3.6 R©, and a Roche lobe for the donor 
of 1.2 Rq. At such a radius, a subgiant expands on a time 
scale longer than R/R ~ 10"'^''yr. If the donor has already 
reduced its mass by mass transfer to the neutron star, 
it Roche lobe will be somewhat smaller (e.g. 0.9 Rq at 
a mass of 0.4 M0), and the expansion time scale longer. 
For conservative mass transfer, in which no mass or an- 
gular momentum is lost from the binary, the time scales 
for the mass transfer and for changes in the orbital pe- 
riod are set by the expansion time scale of the donor, i.e. 
M/M ~ P/P ~ R/R. The time averaged mass transfer 
rate of M ~ 1.4 x 10^^°Mq yr"^ derived above, and the 
limit of the period change P/P < 3 x 10"* yr"^ are both 
compatible with the expansion rate of a subgiant donor 
star. The low mass transfer rate furthermore indicates that 
the time scale for angular momentum loss from the binary 
through magnetic braking is longer than ~ 10^° yr. 

An alternative explanation for an expanded donor star 
would be that the donor still contains a fair amount of 
thermal energy left over from its capture by the neutron 
star (e.g., Verbunt 1994). Since the life time of such a 
source would be set by its thermal time scale, this possi- 
bility is less probable a priori; and the arguments above 
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show that it is not necessary as ah the properties of the 
low-mass X-ray binary may be explained with an ordinary 
subgiant donor. 

7.3. Eclipse profile 

The transition profiles in GRS 1747-312 are remarkably 
stable, when compared to other LMXBs. Ingress and 
egress durations range from -30 to +30% with respect to 
the average. In EXO 0748-676, for instance, the variation 
is 100% (Wolff et al. 2002). The stability may be due to 
either a relative stability of magnetic activity on the sec- 
ondary, or a relative smoothness of the accretion process 
(and associated X-ray irradiation of the secondary). 

We observed a difference in the eclipse duration at 2-5 
keV and 5-20 keV energies and modeled this by assuming 
an isothermal spherical atmosphere in hydrostatic equi- 
librium where the density is an exponential function of 
height. This simple model reproduces the eclipse duration 
difference for a scale height, as projected on the path of 
the neutron star line of sight through the atmosphere, of 
2.0 X 10^ km for the ingress and 2.8 x 10^ km for the 
egress. The model reproduces the curvature on both ends 
of the transition profiles and the spectral change over the 
transition. The true scale height depends on the precise in- 
clination angle. If the primary and secondary are 1.4 and 
0.8 Mq respectively, the inclination angle would be 74?5, 
the secondary radius 1.2 Rq and the true scale height 
1.2 X 10"^ km. This is about 2 times larger than expected, 
but this could be explained by X-ray irradiation of the 
secondary. 

The « 25% difference between the average ingress and 
egress duration has not been reported before in any other 
LMXB cclipscr. Perhaps the relative stability of the tran- 
sitions is what makes this measurement possible in this 
particular source. The difference is indicative of a sim- 
ilar difference in the atmosphere's scale height between 
the leading and trailing hemisphere of the secondary. Day 
et al. (1988) predicted such a scale height difference, in 
the same sense, from a Coriolis force acting on supersonic 
flows on the secondary's surface as invoked by X-ray ir- 
radiation. This will increase the scale height on the trail- 
ing hemisphere while decreasing it on the leading hemi- 
sphere. Whether this model is viable to GRS 1747-312 is 
hard to verify without the detection of the secondary star 
(note that the secondaries have been detected in the well- 
documented eclipsers Her X-1, MXB 1658-298 and EXO 
0748-676). 

7.4. Bursts 

The first-time detection of type-I X-ray bursts miambigu- 
ously proves that the compact object in GRS 1747-312 is a 
neutron star. This leaves AC211 in M15 as the sole bright 
LMXB in a Galactic globular cluster for which no na- 
ture has been determined of the compact object. Finding 
bursts in this object is going to be at least as difficult 



as for GRS 1747-312, because it also is a high-inclination 
system. In fact, the central source is obscured by the ac- 
cretion disk rim all the time (e.g., Ilovaisky et al. 1993) 
and bursts can only be detected indirectly through scat- 
tered X-rays. Another complication for AC211 is that it is 
accompanied by another bright and bursting LMXB at a 
distance of only 2'.'7 (White & Angelini 2001) so that only 
Chandra and XMM-Newton will be able to localize bursts 
accurately enough. 

All seven bursts from GRS 1747-312 are relatively 
short suggesting that the flashes occur in a hydrogen- 
poor/hclium-rich layer which, according to burst theory 
(Fujimoto et al. 1981), must have been formed by sta- 
ble hydrogen fusion. For that to happen, the accretion 
rate must be higher than about 1% of Eddington. This 
value is consistent with the observations. The highest ob- 
served peak flux translates to a luminosity of (2.3j^i f ) x 
Iq38 (^j.g g-1 j^Qj. ^ distance of 9.5^2 5 kpc (Kmilkcrs ct al. 
2003). This value is consistent with the Eddington limit as 
expected for a hydrogen-rich atmosphere. The four bursts 
that were well measured with the PCA have the following 
interesting characteristics: 

1. there is a large range of peak fluxes from 0.2 to 2.1 x 
10"^ erg cm~^s~^; 

2. the two radius expansion bursts have fast rise times 
below 1 s, while the others have slower rise times of at 
least 2 s. In fact, burst no. 1 has a rise and decay of 
similar duration which is atypical for a type-I burst; 

3. one radius-expansion burst (no. 3) has a 1.8 times 
higher unabsorbed bolometric peak flux than the 
other; 

4. the minimum apparent radius of a spherical blackbody 

is as low at 4 km. 

We suspect that at least three of these atypical observa- 
tions are associated with the high inclination angle of the 
binary. This suspicion is motivated by a recent observa- 
tion of MXB 1658-298 by Wijnands et al. (2002). They 
observed two bursts during a period of dipping activity 
which had peak intensities that were four and ten times 
smaller than of most ordinary bursts and which showed 
slow rises that were similar to the decay times. We should 
point out that there is one characteristic which makes the 
two MXB 1658-298 bursts different from the GRS 1747- 
312 ones: they were separated by merely (a record) 50 s 
which suggests that incomplete burning (e.g., Fujimoto et 
al. 1987) may have been important in determining some 
of the burst characteristics. Short recurrence times can 
be excluded for GRS 1747-312: the data coverage around 
the two small bursts without additional bursts is at least 
1200 s. However, it is our impression that the short recur- 
rence time does not explain the slow rise times. No dipping 
was observed around the bursts of GRS 1747-312, but we 
note that the two small bursts were observed at orbital 
phases of 0.87 (burst 1) and 0.49 (burst 4) that are more 
favorable for obscuration effects than for the bright bursts 
(the phase is 0.28 for burst 2 and 0.21 for burst 3). 
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Regarding the factor-of-1.8 peak flux difference be- 
tween the two radius expansion bursts, this is similar as in 
another well-measured system. Galloway et al. (2003) were 
able to gather observations of 61 photospheric radius ex- 
pansion bursts from 4U 1728-34 (=GX 354-0). They found 
a variation in the peak flux by a factor of 1.5 which is cor- 
related with spectral variations of the persistent emission. 
They suggest that a variation in the geometry of the ac- 
cretion disk may be at work, for instance precession of a 
warped disk. 
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